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Abstract 



The Lipkin-Nogami method is generalized to deal with finite range density 
dependent forces. New expressions are derived and realistic calculations with 
the Gogny force are performed for the nuclei 164 Er and 168 Er. The sharp 
phase transition predicted by the mean field approximation is washed out by 
the Lipkin-Nogami approach. A much better agreement with the experimental 
data is reached with the new approach than with the Hartree-Fock-Bogoliubov 
one, specially at high spins. 
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Mean-field theories (BCS, Hartree-Fock and Hartree-Fock-Bogoliubov, HFB) are the 
cornerstone of all microscopic approximations to the nuclear many body problem. The 
success of these approaches is mostly due to their ability to deal with single particle motion 
as well as with the collective motion associated with symmetries Q the latter one by enlarging 
the variational Hilbert space with wave functions that are not eigenstates of the symmetry 
operators. These wave functions are usually constrained to accomplish the symmetries on 
the average; for most symmetries (with the exception of the particle number) this is a 
very satisfactory approach, see Ref. |] for a thorough discussion. In the case of pairing 
correlations, in which we are interested in this letter, the crucial quantities are the number 
of correlated pairs and the level density around the Fermi surface. If these quantities are 
small, and in nuclei they usually are, mean field theories are not enough and one should do 
something better. 

The ideal treatment of pairing correlations in nuclei is particle number projection before 
the variation ||T|. At high spin this theory is rather complicated and up to now it has only 
been applied with separable forces 0. On the other hand, the semi-classic recipe of solving 
the mean field equations with a constraint on the corresponding symmetry operator can 
be derived as the first order result of a full quantum-mechanical expansion (the Kamlah 
expansion) || of the projected quantities. The second order in this expansion takes into 
account the particle number fluctuations and might cure some of the deficiencies of the first 
order approximation. However, full calculations up to second order are rather cumbersome 
and just a simple model calculation has been carried out up to now f|]. Most second order 
calculations have been done using the Lipkin-Nogami (LN) recipe originally proposed in 
Refs. see also the clarifying papers by Quentin et al. || and by Flocard at al. [|J. 

Recently the LN method has been applied to study superdeformed nuclei at high spins 



Continuous symmetries as rotations in any space: coordinate space, gauge space of particle 
number operator, etc, as well as with discrete symmetries, e.g. spatial parity. 
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adding a monopole pairing interaction to the Woods-Saxon plus Strutinsky || and to the 
Skyrme force flIU| . 

Up to now the LN method has been only formulated for the case of a simple separable 
pairing interaction. The purpose of this Letter is to extend such studies to more realistic 
pairing interactions like the one implicit in the finite range and density dependent Gogny 
force where the particle-hole and particle-particle part of the interaction are generated from 
the same force. We have formulated the Lipkin-Nogami method using the Kamlah expansion 
and treated the density dependence consistently. In the new formulation additional terms 
arise in the equations determining the parameters of the theory. 

Let |$) be a product wave function of the Hartree-Fock-Bogoliubov type, i.e. a particle 
number symmetry violating wave function. We can generate an eigenstate of the particle 
number \^n) by the projection technique [[J 

\y N ) = |$) = -1 [** dW^-^+lQ). (1) 



2tt Jo 

The particle number projected energy is given by 

N (* n \H\* n ) = ^dtPe-^^He^^) = h(4>) 

proj (^nI^n) " J 27r (i(/)e-^($|e^|$) J 2n d(f)e-^ N n(4>) 1 ' 

where we have introduced the hamiltonian-, h(<p) = ($|ife 1 ^ |3>), and norm-, n(0) = 
($|e* <?iiV | ( l ) ), overlaps. In the case of large particle numbers and strong deformations in 
the gauge space associated to N, one expects the h(<p) and n(<p) overlaps to be peaked at 
= and to be very small elsewhere in such a way that the quotient h(4>) /n((f>) behaves 
smoothly. One can make an expansion of h(4>) in terms of n(4>) in the following way |J 

M 

h{<f>) = E h m M m n(<P) (3) 

m=0 

where we have introduced the Kamlah operator J\f = \-§g — ($|iV|$) which is a representation 
of the particle number operator in the space of the parameter 0. The expansion coefficients 
h m are determined by applying the operators l,Jx, ...,J\f M on Eq. (^) and taking the limit 
cf) —>■ 0. From now on we shall use the shorthand notation (A) = and AiV = N—(N). 



The projected energy to second order is 



E™ oj = (H) h 2 ((AN) 2 ) + h(N-(N)) + h 2 (N-(N)) 2 , 



(4) 



with h\ and h 2 given by 

(HAN) - h 2 ({AN) 3 ) 



h 2 



((ANf) 

((H - (H))(ANf) - (HAN) ((AN) 3 ) /((AN) 2 ) 



(5) 
(6) 



((AN) 4 ) - ((AN) 2 ) 2 - ((AN) 3 ) 2 1 ((AN) 2 ) 
In a full variation after projection method one should vary Eq. (fj). In the Lipkin-Nogami 
prescription, however, the coefficient h 2 is held constant during the variation; the resulting 
equation is much simpler, one gets 



(7) 



with hi determined by the constraint 



(N) = N. 



(8) 



If there are additional constraints, for example the angular momentum, one just has to 
substitute H by H' = H — uJ x in all equations above and to add the constraint (J x ) = 
[1(1 + 1) - (J z 2 )} 1/2 to the one of Eq. (§). 

Now we would like to generalize the formulae above to density dependent forces, like the 



Gogny force [11, Eg] which has a term proportional to p a ( ri % r2 , 



In this case the density termQ causes a dependence on <fi of the hamiltonian [T3|]. The 
Kamlah expansion (H) provides in this case the following equation system 



1 d 2 H 



+ 2 AN 



<t>=Q 



lldH 

\ i d(f> 

IdH 

i deb 

<t>=o 



0=0 



+ HAN^j = h((AN) 2 ) + h 2 ((AN) 3 ) 



(9) 



+ H(AN) 2 ) = (H)((AN) 2 ) 



+ h((AN) 3 ) + h 2 (((AiV) 4 ) - ((AiV) 2 ) 2 



2 The dependence of the density with 4> is given by p(r) = (<&|c^(r)c(r)e^ Ar |$)/($|e i ^ >Ar |<I>). 
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which determine the coefficients hi and hi- 

From now on we proceed as in the non-density dependent case, i.e. we have to solve 
Eq. ( [7| ) with the constraint (§) but with the coefficient hi and hi provided by the equation 
system flj^). 

We have applied this formalism to study high spin states with the Gogny force in two rare 
earth nuclei: 164 Er as an example of strong back-bender and 168 Er as a non-back-bender. We 
shall refer to this calculations as cranked-HFB-LN (CHFBLN). In all our results we shall also 
present the ones obtained with the plain cranked HFB theory (CHFB) Jl4|]. We are using 
the standard DS1 parametrization set fl2|l . Of course, one could ask if the parametrization 
of the force should be changed, for it was adjusted for plain mean field calculation. Since 
this is the first investigation in this direction we shall use the standard parametrization, 
further investigations will be reported in the future. 

In Fig. |l|a we show the pairing energy in the CHFB approach and in the CHFBLN 
approach versus the angular momentum. For 164 Er in the CHFB approach we observe for 
the neutron system, first, a strong Coriolis antipairing effect which diminishes the neutron 
pairing correlations, later on the crossing of the ground state band with a two neutron 
aligned band -see below- causes the quenching of the pairing correlations at / ~ 18ft. The 
proton system, on the other hand, behaves very smoothly until / = 18ft. From I = 18h 
up to / = 28ft we observe the typical Coriolis antipairing effect reduction, which is not as 
strong as for the neutron system because the intruder orbital in this case is a ftn at variance 
with the ii3 of the neutrons. In the CHFBLN results the same Coriolis antipairing effects 
are observed but no superfluid to normal fluid phase transition is found. We also realize 
that the LN term has a larger effect on the proton system than on the neutron one. For 
168 Er, Fig. |I|b, again the neutron phase transition is washed out in the CHFBLN approach 
and a larger increase in the pairing energies of the proton system than in the neutron one 
is obtained; this may have to do with the different intruders for both systems. 

The most relevant deformation parameters are f3 and 7 (we define them as in Ref. [|H)|]). 
Their angular momentum dependence is displayed in Figs. |I|c and |l]d. For 164 Er, in the 
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CHFB approach we first observe a rather constant value of the deformation parameters (3 
until / = 12h, from this point on and until / = 28h we find a decrease in f3. This anti- 
stretching effect is caused by the Coriolis force. The CHFBLN approach differs from the 
CHFB one in the spin range / = 10h till I = 18h, where the neutron pairing collapse take 
place. Along the aligned band the /3 values are again very similar in both approaches. In the 
CHFB approach, at I = Oh the nucleus 164 Er is axially symmetric (7 = 0), then 7 increases 
up to 8 degrees at / = 18h, later on it decreases very slowly. In the CHFBLN approach, the 
nucleus remains axially symmetric (7 = 0) at all spin values. In the nucleus 168 Er we observe 
a similar behavior when we compare the /5-values of the CHFBLN and the CHFB at the 
spin range I = 12 — 18h. From 1 = 22 — 28% we find a larger decrease in the /3-values of the 
CHFBLN as compared with the CHFB. The reason for this behavior can be found in Fig. ^d; 
due to the smaller value of the moment of inertia in the CHFBLN approach as compared 
with the CHFB one, larger values of the cranking frequency are needed to produce the same 
angular momentum. These larger values of the cranking frequency causing an stronger anti- 
stretching effect on the /3-values of the CHFBLN approach. The 7-values of 168 Er are close 
to zero in both approaches. 

The E2 transition probabilities and the gyromagnetic factors have been calculated in 
the cranking approximation |14[]. Since our configuration space is large enough (11 oscillator 



shells) no effective charges have been used in the calculations. In Fig. |2|a we show the 
reduced transition probabilities along the Yrast band versus the angular momentum and 
the experimental ones for the nucleus 164 Er. At spin values up to / = lO/i our theoretical 
results are in good agreement with the experimental data. For spin values 12, 14 and 
16fr, corresponding to the band crossing, we are not able to reproduce the zig-zag behavior 
of the experimental data. This result is not surprising since we know that the cranking 
approximation is not good in the band crossing. The decrease of the CHFBLN results 
as compared with the CHFB ones is due to the fact that in the CHFBLN approach the 
nucleus remain axially symmetric at this spin values while in the CHFB approach it does 
not. Concerning 168 Er, Fig. |2|b, neither CHFB nor CHFBLN are able to reproduce the zig- 
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zag at medium spins (this behavior is not understood, to our knowledge, in any theory). 
The smaller values of the CHFBLN at high spins are due to the smaller /3-values of this 
approach. 

To investigate the alignment processes in these nuclei we shall study the gyromagnetic 
factors. In the calculation the free values of the orbital and spin gyromagnetic factors have 
been used and no rotor contribution has been considered. In Fig. ||]c we display the 
theoretical g— factors g,g p and g n as well as the experimental g— factor at / = 2% [|l^] for 
164 Er. From the pattern of g we can conclude that at low spins we have a smooth neutron 
alignment, at medium spins and up to I = 18h we have strong neutron alignment and for 
higher spins we observe proton alignment. The tendency is qualitatively the same in both 
approximations, the CHFBLN displaying a sharper behavior though. For 168 Er, Fig. 0d, we 
obtain in the CHFB (CHFBLN) an smooth neutron alignment up to spin 16h (24h), later 
on proton alignment. The agreement with the known experimental data (/ = 6, 8 and 10h) 
is excellent. 

In Fig. |^a we display transition energies versus the angular momentum. For 164 Er the 
agreement with the experimental data is good at low spins in the CHFB approach while in 
the band crossing region we see that the crossing is not as sharp as in the experiment -this is 
a well known drawback of the cranking approximation. In the high spin part we get smaller 
values for the transition energies than in the experiment. In the CHFBLN approach at low 
and medium spins the agreement with the experiment is better than in the CHFB approach. 
In the backbending region the results are not good again, but in the very high spin region 
the agreement with the experiment is very good at variance with the CHFB approximation. 
Concerning 168 Er, Fig. ^jp, the CHFB results are very low as compared with the experiment; 
the CHFBLN ones, however, are in excellent agreement with the experiment. 

In Fig. |3]c, d we display the moments of inertia, versus the square of the angular frequency. 
In 164 Er in the CHFBLN approach we obtain at low spins a smaller value than in the CHFB 
due to the larger pairing correlations in better agreement with the experiment. In the band 
crossing we obtain a clear back-bending although shifted in a few units as compared with 
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the experiment indicating, perhaps, that angular momentum projection is important in this 
region. For spin values I = 18, 20, 22h, i.e. on the aligned band, the agreement with the 
experiment is excellent [], at variance with the simpler CHFB approach. Concerning 168 Er, 
Fig. |3]d, the results of the CHFBLN approach are much better than the ones of CHFB, we 
are at low spins somewhat lower than in the experiment, but at medium and high spins the 
agreement is very good. We would like to stress that no J C) the core moment of inertia, has 
been assumed. 

In conclusion, for the first time, we have formulated the Lipkin-Nogami approximation 
for density dependent forces. We have performed numerical calculations with the finite range 
density dependent Gogny force for two nuclei, the theoretical results with this approximation 
are in much better agreement with the experiment than the plain HFB calculations. 

This work was supported in part by DGICyT, Spain under Project PB9 1-0006. One 
of us (A.V.) would like to thank the Spanish Ministerio de Asuntos Exteriores for financial 
support through an ICI grant. 



3 Notice that the experimental results end at spin / = 22h {I = 16h) for 164 Er ( 168 Er) while the 
theoretical ones go up to / = 26h. 
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FIGURES 

FIG. 1. Upper panels: Pairing energies versus the angular momentum. Proton (neutron) 
values are represented by triangles (inverted triangles). Lower panels: Deformation parameters (5 
(circles, scale on the left axis) and 7 (squares, scale on the right axis) versus angular momentum. 
Full (empty) symbols correspond to CHFB (CHFBLN). 

FIG. 2. Upper panels: Reduced transition probabilities B(E2) (in units of (eb) 2 ) along the Yrast 
band as a function of the angular momentum. Full (open) circles stand for the CHFB (CHFBLN) 
results and full squares for the experimental data ( for 164 Er and 
panels: Gyromagnetic factors (in units of //jv of the Yrast states versus the angular momentum; g 
is represented by diamonds, g p by circles and g n by triangles. Full (open) symbols correspond to 
CHFB (CHFBLN). The experimental data ( @ for 164 Er and @ for 168 Er) are represented by 
full squares. 

FIG. 3. Upper panels: the gamma-ray energy AEj = E{I) — E(I — 2) as a function of the 
angular momentum. Full (open) circles stand for the CHFB (CHFBLN) results and full squares 
for the experimental data ( [l(| for 164 Er and ^L7§ for 168 Er). Lower panels: the moment of inertia 
J = {21 — 1)/ AEj versus the square of the angular frequency. The meaning of the symbols is the 
same as in the upper panels. 
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